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Abstract 

Shortly  after  initiating  this  program,  we  focused  our  effort  to  develop  thermodynamic 
descriptions  of  ternary  Ni-Cr-Ru  and  Ni-Al-Ru  including  the  constituent  binaries  when 
needed  using  the  traditional  Calphad  approach.  In  addition,  we  extended  our  effort  to 
include  the  use  of  the  Cluster/Site  Approximation  (CSA)  to  describe  the  fee  phases  in  the 
disordered  and  ordered  states  such  as  the  prototype  Cu-Ag-Au  in  2003/2004  and  then  y/y 
in  real  ternary  systems  in  2004/2005.  In  the  meantime,  we  began  to  explore  the 
possibility  of  using  CSA  to  calculate  interphase  boundary  (IPB)  energies  of  coherent 
interfaces  such  as  those  for  y/y  in  prototype  Cu-Au  binary  and  then  real  binaries  such  as 
Ni-Al.  The  preliminary  results  on  CSA-calculated  prototype  Cu-Ag-Au  diagrams  and 
IPB  energies  on  prototype  Cu-Au  and  Ni-Al  were  presented  at  the  2004  annual  review 
meeting  at  Wintergreen,  VA.  We  subsequently  discussed  with  Dr.  C.  Hartley  about  our 
IPB  energy  effort  and  were  encouraged  to  continue  this  research  effort  as  part  of  our 
program. 

Developing  a  thermodynamic  description  of  Ni-Al-Ru  has  been  truly  challenging  using 
the  limited  and  often-inconsistent  experimental  data  reported  in  the  literature.  Obtaining 
a  thermodynamic  description  for  an  alloy  system  is  not  purely  a  mathematical  fit  of  the 
experimental  data  to  a  set  of  mathematical  equations.  Rather  the  obtained  parameters 
must  be  physically  reasonable  in  terms  of  the  enthalpy  and  entropy  quantities.  We  had 
every  reason  to  believe  that  the  phase  equilibrium  data  reported  in  the  literature  could  be 
in  error.  Nevertheless  we  tried  our  best  to  obtain  a  set  of  thermodynamic  parameters  for 
all  the  phases  in  this  alloy  system.  Based  on  the  calculated  phase  diagrams  using  the 
preliminarily  obtained  thermodynamic  description,  we  identified  a  limited  number  of 
alloy  compositions  for  investigation  either  to  verify  the  calculation  or  more  likely  in  this 
case  to  allow  us  to  obtain  an  improved  description!  The  experimental  results  indicate  that 
this  indeed  is  the  case.  An  improved  thermodynamic  description  of  this  technologically 
important  ternary  was  then  obtained  both  in  the  traditional  Calphad  approach  and  using 
the  CSA  to  describe  all  the  fee  related  phases. 

We  successfully  extended  the  CSA  to  calculate  phase  diagrams  beyond  binaries  first  to  a 
prototype  ternary,  next  to  Ni-Al-Cr,  and  lastly  quaternary  Ni-Al-Cr-Re.  The  CSA 
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calculated  phase  diagrams  are  in  accord  with  available  data.  Moreover,  few  parameters 
are  needed  when  compared  with  previous  descriptions  using  the  point  approximation.  In 
addition,  the  CSA-calculated  metastable  fee  phase  diagrams  and  the  fee  +  liquid  phase 
diagrams  are  what  one  expects.  We  also  demonstrated  that  the  CSA  possesses 
comparable  accuracy  when  compared  with  CVM  in  phase  diagram  calculation  but  is 
computationally  much  less  demanding. 

The  advantage  of  computational  speed  makes  CSA  a  viable  model  to  calculate  not  only 
multicomponent  phase  diagrams  but  also  coherent  interphase  boundary  (IPB)  energies. 
We  first  demonstrated  its  suitability  for  calculating  coherent  IPB’s  for  prototype  ordering 
and  phase-separating  alloy  systems  using  the  CVM-calculated  values  as  experimental 
data.  We  next  extended  this  approach  to  calculate  IPB  energies  between  y-(Ni)  and  f- 
(NhAl)  in  the  Ni-Al  and  (Al)  and  (AI3U)  in  Al-Li,  in  accord  with  the  CVM-calculated 
values  and  available  experimental  data  within  the  uncertainties  of  the  measured  values. 
In  addition,  we  also  extended  the  IPB  energy  calculation  to  anti-phase  boundary  (APB) 
energies  in  prototype  binaries  using  a  similar  approach  as  has  been  done  for  IPB  energies. 
These  results  indicate  that  one  can  apply  CSA  to  calculate  the  coherent  y/y’  IPB  and  APB 
energies  in  multicomponent  Ni-based  superalloys. 

Introduction,  Research  Objective  and  Progress 

Improvement  of  future  generation  turbine  engine  propulsion  systems  requires 
development  of  better  materials  that  are  more  durable  under  the  harsh  conditions  of  high 
temperature  and  cyclic  stress.  Traditional  development  cycles  for  turbine  engine  Ni- 
based  superalloys  can  last  more  than  10  years.  This  is  partially  due  to  the  difficulty  of 
finding  an  alloy  chemistry  that  provides  the  desired  mechanical  properties  and  material 
stability.  The  best  way  to  accelerate  the  development  of  new  materials  is  to  develop 
advanced  computational  thermodynamics  as  an  essential  tool  to  provide  a  blueprint  for 
alloy  development  and  process  optimization.  Computational  thermodynamics  enables  us 
to  calculate  multicomponent  phase  diagrams  that  provide  the  road  maps  for  selecting 
alloy  chemistry  and  optimizing  processing  parameters  for  solidification  and  heat 
treatments.  The  objective  of  the  present  study  is  to  adopt  the  Calphad  approach  to 
develop  thermodynamic  descriptions  of  Ni  alloys  containing  elements  such  as  Al,  Cr,  Ru, 
etc  and  to  simultaneously  develop  CSA  as  an  enabling  model  to  describe  the 
thermodynamics  of  the  fee  phases  such  as  y  and  y’  in  Ni-based  superalloys. 

Phase  Diagram  calculation  and  experimental  investigation  of  the  Ni-Al-Cr-Ru  System 

The  phase  equilibria  in  the  Ni-Al-Ru  system  have  been  studied  by  a  number  of 
investigators  [80Tsu,  85Cha,  85Pet,  86Cha,  97Har,  97Hor,  97Horl,  98Hor,  OOHoh]. 
Most  of  these  investigations  have  focused  on  solid-state  phase  equilibria.  The  primary 
phases  of  solidification  reported  by  Hohls  et  al.  [OOHoh]  and  Homer  et  al.  [97Horl] 
focused  on  the  Al-rich  corner  whereas  the  interest  here  is  primarily  on  the  Ni-rich  comer. 
Obtaining  a  thermodynamic  description  of  this  ternary  has  been  extremely  challenging 
due  to  the  complex  phase  equilibria  and  a  lack  of  extensive  and  internally  consistent 
experimental  data.  For  the  phase  equilibrium  data,  we  have  relied  primarily  on  the 
isothermal  sections  obtained  by  Chakravorty  and  West  [86Cha]  and  we  did  try  to  use 
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other  earlier  published  data  in  our  effort  to  obtain  a  thermodynamic  description  but  it  was 
impossible  to  achieve  internal  consistency.  The  only  enthalpy  of  formation  data  available 
in  this  system  is  that  of  Jung  and  Kleppa  for  (3-RuAl  [92Jun].  This  value  was  adjusted  to 
be  more  negative  since  the  enthalpy  of  formation  for  (J-NiAl  reported  by  Nash  and 
Kleppa  [OINas]  was  less  negative  than  the  accepted  data  reported  by  Henig  and  Lukas 
[75Hen],  We  first  used  the  traditional  Calphad  approach  to  model  this  system,  i.e.,  the 
disordered  y-(Ni)  phase  and  the  ordered  Y-(Ni,Ru)3Al  phase  were  treated  as  two  different 
phases.  During  the  optimization  of  model  parameters,  we  found  that  it  was  very  difficult 
to  maintain  the  phase  relationships  at  1273K  reported  by  Chakravorty  and  West  [86Cha], 
The  y-(Ni,Ru)3Al  phase  tends  to  be  in  equilibrium  with  the  5-(Ru)  phase.  In  order  to 
achieve  equilibrium  between  P  and  y-(Ni),  we  needed  to  use  a  very  large  negative  ternary 
interaction  parameter  for  the  y-(Ni)  phase.  Even  worse  situation  was  met  when  we  tried 
to  model  this  system  using  the  more  physical  realistic  CSA  model  to  describe  the  two 
fcc-based  phases  (y  and  y1)  simultaneously.  In  this  case,  the  Y-(Ni,Ru)3Al  phase  was 
always  in  equilibrium  with  the  5-(Ru)  phase  at  1273K  no  matter  how  we  changed  the 
model  parameters.  We  had  every  reason  to  believe  that  the  phase  equilibrium  data 
reported  in  the  literature  are  likely  in  error. 

Experiments  were  carried  out  in  order  to 
confirm  the  phase  equilibria  at  1273K  in 
the  Ni-rich  corner  as  well  as  the  primary 
phases  of  solidification.  Based  on  the 
calculated  phase  diagrams  using  the 
preliminarily  obtained  thermodynamic 
description,  we  identified  a  number  of 
alloy  compositions  for  investigation. 

The  alloy  samples  were  prepared  by  arc¬ 
melting  using  high  purity  metals  under  a 
partial  pressure  of  argon  and  then  cut 
into  portions.  One  portion  was  for 
investigation  in  the  as-solidified 
condition.  The  remaining  portion  was 
sealed  in  quartz  tube  under  a  partial 
pressure  of  argon  for  annealing  at  1273K  for  30  days  followed  by  quenching  in  water. 
The  as-cast  and  annealed  samples  were  characterized  by  XRD,  SEM  and  EPMA.  Several 
annealed  samples,  in  the  y  +  o  +  P2  three-phase  region  according  to  Chakravorty  and 
West,  show  the  existence  of  y'  +  a  +  P2,  consistent  with  the  CSA  calculation  but  in 
contradiction  with  the  data  of  Chakravorty  and  West.  Fig  1  shows  the  XRD  patterns  for 
the  alloy  Ni-20Ru-27Al  (at.%)  after  30-day  annealing  at  1273K.  These  data  clearly 
indicate  that  Y-(Ni,Ru)3Al  and  8-(Ru)  are  in  equilibrium  with  each  other  rather  than  P  and 
y-(Ni).  In  Chakravorty  and  West’s  experimental  work,  the  alloys  were  annealed  for  only 
10  days  at  1273K  after  10-day  annealing  at  1523K.  At  1523K  there  are  three  stable 
phases  in  the  sample:  y  +  p2  +  a.  When  the  sample  was  cooled  down  to  1273K,  the  y 
phase  did  not  have  sufficient  time  to  transform  to  Y  due  to  the  small  driving  force  from 
the  disordered  y  phase  to  ordered  Y  phase.  The  results  we  obtained  when  the  samples 
were  annealed  for  one  week  at  1273K  were  the  same  as  those  obtained  by  Chakravorty 


after  30-day  annealing  at  1273K 
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and  West,  confirming  that  a  longer  annealing  time  is  necessary  to  achieve  equilibrium! 


Fig.  2  Comparisons  between  the  CSA-calculated  Ni-Al-Ru  isothermal  sections  with  experimental 
data  (a)  1523  K  (1250°C),  and  (b)  1273  K  (1000°C) 


To  make  a  comparison  between  the  calculated  results  and  the  available  experimental  data 
[86Cha],  two  isothermal  sections  were  calculated  at  1523  and  1273K.  As  shown  in  Fig. 
2(a),  the  calculated  three-phase  equilibria  of  (3  +  y-(Ni)  +  8-(Ru)  and  (3  +  y-(Ni)  +  y1- 
(Ni,Ru)3Al  are  in  reasonable  accord  with  the  measured  compositions  of  the  co-existing 
phases.  In  the  isotherm  at  1273K  shown  in  Fig.  2(b),  the  agreement  between  calculated 


Fig  3.  Model-calculated  liquidus  projection 
of  Ni-AI  -Ru  ternary  system 


Fig  4.  Comparison  of  calculated  Ni-Cr-Ru  1250 
°C  isotherm  with  experimental  data  [85Cha] 


and  our  own  experimental  results  is  comparable.  As  shown  in  Fig.  3,  the  calculated 
liquidus  projection  shows  a  large  composition  region  where  the  primary  phase  of 
solidification  is  (3,  due  to  its  rather  high  melting  point.  Below  50  at%  Al,  there  is  a  saddle 
point  in  the  calculated  liquidus,  i.e.  a  critical  tie  line  of  L  +  a  +  y  at  this  temperature.  In 
addition  there  exist  two  type  II  four-phase  invariant  equilibria:  IF  (L  +  y1  =  |3  +  y)  at 
1592K  and  4  at%  Ru,  and  II2  (L  +  5  =  [3  +  y)  at  1575K  and  61  at%  Ni,  21  at%  Al, 
respectively.  Above  50  at%  Al,  there  is  a  cascade  of  invariant  equilibria  towards  the  Al- 
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rich  comer  with  decreases  in  temperature. 

The  phase  equilibria  in  Ni-Cr-Ru  ternary  are  much  simpler.  Figure  4  shows  a  model- 
calculated  1250°C  isotherm  of  Ni-Cr-Ru.  It  is  in  agreement  with  the  data  of  [85Cha]. 
The  Al-Cr-Ru  ternary  is  far  away  from  the  Ni  comer,  so  a  thermodynamic  description  of 
this  ternary  was  obtained  by  extrapolation  from  the  three  constituent  binaries. 

The  thermodynamic  description  of  the  Ni-Al-Cr-Ru  system  was  obtained  from  those  of 
the  constituent  binary  and  ternaries  by  extrapolation.  As  shown  in  Fig.  5,  the  calculated 
isothermal  sections  with  75  at.  %  of  Ni  are  in  good  agreement  with  the  experimental  data 
[85Chal].  This  quaternary  description  takes  into  account  the  experimental  data  which 
are  judged  to  be  reliable  and  is  believed  to  offer  an  understanding  of  the  phase  stability  in 
the  Ni-rich  comer  and  a  useful  guide  for  further  experimental  investigation  in  other 
compositional  regions. 


Fig  5.  CSA-calcuIated  isothermal  sections  with  75  at.%  of  Ni  with  experimental  data[85Chal] 

(a)  1273K  and  (b)  1523K 


Phase  diagram  Calculation  -  Cluster/Site  Approximation 

The  Cluster/Site  Approximation  (CSA)  takes  into  account  short-range  order  (SRO), 
which  is  essential  to  satisfactorily  describe  the  thermodynamics  of  order/disorder 
transitions  such  as  the  transition  from  y -(fee)  to  y'-{L\2)  in  multicomponent  Ni-based 
superalloys.  LI 2  is  an  ordered  form  of  the  fee  structure.  As  noted  earlier  CSA  offers 
computational  advantages  over  the  Cluster  Variation  Method  (CVM)  while  retains 
comparable  accuracy  when  compared  with  CVM-calculated  phase  diagrams.  This  makes 
the  CSA  a  practical  model  for  calculating  phase  diagrams  of  real  multicomponent  alloy 
systems.  We  first  demonstrated  the  suitability  of  CSA  in  calculating  prototype  binary 
phase  diagram  using  CVM-calculated  diagrams  as  the  reference  [990at,  OlZha,  04Cha] 
and  then  extended  this  approach  to  calculate  the  technologically  important  binary  Ni-Al 
phase  diagram  [03Zha,  04Cha], 
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In  addition,  we  have  successfully  extended  this  approach  to  calculate  the  phase  diagrams 
of  the  prototype  coherent  ternary  Cu-Ag-Au  system  [2005Cao],  again  using  the  CVM- 
calculated  diagrams  as  the  reference  [80Kik],  Good  agreement  was  again  achieved. 
However,  it  is  noteworthy  to  point  out  that  the  CSA-calculated  isotherm  at  240°C  as 
shown  in  Fig.  6(b)  differs  from  that  calculated  with  CVM  shown  in  fig.  6(a).  Figure  6  (b) 
does  not  show  the  existence  of  the  Ll2  phase  in  the  Au-rich  corner.  Prof.  C.  Colinet  of 
LTPCM-ENSEEG,  Saint  Martin  d’Heres,  France  subsequently  re-calculated  this  isotherm 
using  the  CVM  and  indeed  found  the  existence  of  the  LI 2  phase  in  the  Au-rich  corner.  In 
the  original  CVM  calculation,  it  was  necessary  to  know  the  existence  of  this  Ll2  phase  in 
the  Au-rich  comer  so  as  to  estimate  initial  values  to  calculate  the  compositional  stability 
of  this  phase.  On  the  other  hand,  we  use  the  Pandat  software,  which  is  able  to 
automatically  find  the  lowest  Gibbs  energy  and  thus  obtain  the  most  stable  phase 
automatically  [2003Che,  2002Che,  2001  Che]. 


Fig  6.  Comparison  of  prototype  Cu-Ag-Au  phase  diagrams  calculated  by  CVM  (a)  and  CSA  (b)  at  240  °C 


Subsequent  to  successfully  extending  the  CSA  from  binaries  to  calculate  prototype 
coherent  phase  diagrams  of  Cu-Ag-Au,  we  applied  this  approximation  to  the  fee  phases 
in  the  Ni-Al-Cr  ternary  system.  Since  the  fee  phases  in  Ni-Al  binary  had  already  been 
modeled  using  the  CSA  by  Zhang  et  al.  [03Zha],  we  adopted  their  description.  The 
model  parameters  for  the  fee  phases  in  the  other  two  constituent  binaries  were  obtained 
from  the  experimental  data  assessed  in  the  literature  [05Caol,  99Hua],  Figures  7(a)  and 
(b)  show  comparisons  between  the  CSA-calculated  isotherm  of  Ni-Al-Cr  at  1273K  and 
the  CSA-calculated  y-solvus  curves  as  a  function  of  temperature  with  experimental  data. 
These  y-solvus  curves  were  calculated  at  constant  values  of  the  mole  fraction  of  Al 
varying  from  0.09  to  0.17.  The  calculated  phase  boundaries  as  shown  in  Fig.  7(a)  are 
consistent  with  the  experimental  data  of  Ochiai  et  al.  [830ch]  and  the  calculated  solvus 
of  the  y-phase  shown  in  Fig.  7(b)  is  also  in  accord  with  the  data  of  Hong  et  al.  [89Hon], 
Many  other  comparisons  between  calculation  and  experimentation  were  made  and  the 
agreement  is  similar  to  those  shown  in  Fig.  7(a)  and  (b)  but  not  presented  here. 
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In  addition  to  be  able  to  calculate  the  stable  binary  and  ternary  phase  diagrams,  it  is 
interesting  to  use  CSA  to  calculate  metastable  phase  diagrams  such  as  the  phase  diagrams 
involving  only  the  fee  phases  and  only  the  fee  and  liquid  phases.  For  instance  the 
topology  of  the  CSA-calculated  metastable  binary  Ni-Al  phase  diagram  for  the  fcc-based 
phase  is  the  same  as  that  calculated  using  the  “first  principle”  combined  with  CVM 
[88Car,  92Pas]  to  be  presented  later.  On  the  other  hand,  Shockley  had  shown  that  the 
calculated  phase  diagram  involving  only  fcc-based  phases  based  on  the  point 
approximation  does  not  have  the  correct  topology  [38Sho],  We  now  present  the  CSA 
model,  used  to  calculate  the  stable  1273  K  isotherm  shown  in  Fig:  7(a),  to  calculate  the 
1273K  isotherm  involving  only  the  fee  phases,  as  shown  in  Fig.  8(a).  As  shown  in  this 
diagram,  the  calculated  y  +  y’  (or  A1  +  LI2)  phase  equilibria  in  the  Ni-rich  comer  are 
nearly  the  same  as  anticipated  since  they  correspond  to  the  stable  equilibria  as  can  be 
seen  from  Fig.  7(a).  Moreover,  the  calculated  metastable  phase  equilibria  away  from  the 
stable  region,  although  different  in  these  two  diagrams,  are  consistent  with  the  phase  rule. 
In  other  words,  the  extensions  of  the  phase  boundaries  to  the  metastable  region  show  the 


Ni 


Figure  7.  Comparison  of  calculation  with  experimental  data  of  Ni-AI-Cr  ternary  system 
(a)  1273K  (1000°C)  isothermal  section  and  (b)  the  y-solvus  at  constant  values  of  A1 


Mole  fraction  Cr  Mole  fraction  Cr 


Figure  8.  Calculated  metastable  isotherms  at  1273K  (1000°C):  (a)  Phase  equilibria  involving 
only  the  fee  phases  and  (b)  Phase  equilibria  involving  only  the  fee  and  the  liquid  phases. 
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existence  of  two-phase  fields  of  LI2  +  Llo,  A1  +  Llo  as  well  as  a  miscibility  of  the  A1 
phase  with  the  corresponding  three-phase  field  toward  the  Al-rich  corner.  For  Fig.  8(b), 
the  liquid  phase  appears  in  the  Al-rich  comer  with  corresponding  two-  and  three-phase 
fields  involving  the  liquid  phase. 

It  is  worth  noting  that  Zhang  et  al.  [2003Zha]  had  shown  a  miscibility  gap  occurring  in 
the  Al  phase  in  the  Al-rich  side  of  the  calculated  metastable  fee  Al-Ni  binary  diagram  as 
shown  in  Fig.  9(a).  To  the  best  of  our  knowledge  none  of  the  descriptions  of  Ni-Al 
obtained  from  the  Calphad  approach  so  far  is  able  to  do  so  [2003Zha,  2004Cha],  We 
present  this  binary  metastable  phase  diagram  of  Ni-Al  at  this  point  because  we  wish  to 
present  an  isopleth  of  Ni-Al-Cr  at  a  constant  value  of  2  at%  Cr  shown  in  Fig.  9(b) 
[2001Caol]  and  make  a  comparison  of  these  two  diagrams.  Figure  9(a)  shows  that  the 


Figure  9  CSA-calculated  metastable  fcc(Al)  phase  diagrams:  (a)  Binary  Ni-Al  [2003Zha] 
and  (b)  An  isopleth  of  Ni-Al-Cr  with  constant  value  of  2  at  %  Cr  [2005Caol] 


Al(fcc)  phase  is  stable  at  high  temperatures.  With  decreasing  temperature,  the  Llo 
appears  at  50  at%  Ni,  the  LI2  phase  at  75  at%  Ni,  and  two  Al  (fee)  phases  toward  the  Al- 
rich  side.  The  calculated  isopleth  of  Ni-Al-Cr  with  2  at%  Cr  has  similar  features  except 
(i)  a  three-phase  field  of  Al  +  Al  +  Llo  appears  between  the  2  two-phase  fields  of  Al  + 
Al  and  Al  +  Llo  and  (ii)  both  the  Llo  and  LI 2  phases  no  longer  melt  congruently  as 
should  be  the  case.  These  results  are  expected  since  a  binary  invariant  such  as  the  binary 
monotectoid,  Al  +  Al  +  Llo,  becomes  a  tie-triangle  in  the  ternary  region  over  a  range  of 
temperatures.  Of  course  there  is  also  another  three-phase  equilibrium  of  Al  +  Llo  +  LI2 
also  shown  in  this  diagram  at  ~  1470K  and  ~  60  at%  Ni.  Many  other  diagrams  could  be 
calculated  such  as  the  liquidus  projection  of  the  metastable  phase  diagram  involving  only 
the  fee  (Al)  and  liquid  phases.  Additional  metastable  diagrams  were  presented 
elsewhere  [05Caol]  and  are  not  repeated  here. 

We  have  successfully  using  CSA  to  calculate  the  quaternary  Ni-Al-Cr-Re  system.  The 
preliminary  results  presented  in  Figs.  10(a)  and  (b)  show  excellent  agreement  between 
calculation  and  experimentation  is  achieved. 
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Fig  10.  The  CSA-calculated  compositions  of  (a)  y  and  (b)  y'  of  the  Ni-Al-Cr-Re  alloys  with  1 .5  at%  Re 

and  xNj-xAi  =0.625  at  1313K. 


Interphase  Boundary  (IPB)  Energy  Calculation  using  CSA 

In  view  of  the  relative  computational  simplicity  of  the  CSA  and  the  fact  that  it  is  able  to 
calculate  phase  diagrams  of  binary  and  higher  order  alloys  system,  we  have  been  first 
exploring  the  use  of  the  CSA  to  calculate  coherent  IPB  energies  for  the  fcc/Ll2  phases. 
Our  motivation  was:  if  this  approach  were  shown  to  be  a  suitable  one,  we  could  then 
extend  the  use  of  this  approach  to  calculate  IPB  energies,  for  instance,  between  y  and  y’ 
first  in  binary  Ni-Al  and  then  in  multicomponent  Ni-based  superalloys.  These  quantities 
are  urgently  needed  for  describing  micro  structure  evolution  when  an  alloy  undergoes 
thermal  treatment  using  kinetic  modeling  such  as  the  phase-field  approach.  We  first  used 
the  CSA  to  calculate  the  IPB  energies  of  binary  immiscible  alloys  and  then  those  between 
the  prototype  fcc-(Cu,  Au)  and  the  ordered  Ll2-(Cu3Au)  phases.  The  calculated  results 
for  the  immiscible  binary  case  agree  with  those  obtained  by  Kikuch  and  Chen  [95Kik] 
using  the  CVM  but  deviate  from  those  obtained  by  Lee  and  Aaronson  [80Lee]  using  a 
regular  solution  model.  The  results  have  been  presented  elsewhere  and  repeated  here 
[2005Cao2], 

Comparisons  of  the  CSA-calculated  coherent  IPB  energies  between  the  fcc-(Cu,  Au)  and 
the  L12-(Cu3Au)  phases  with  the  CVM-calculated  values  [79Kik]  are  presented  in  Figs. 
11(a)  and  (b).  In  both  cases  the  tetrahedron  approximation  was  used.  Our  calculation 
was  carried  out  following  the  approach  of  Kikuchi  and  Cahn  but  using  the  CSA  instead 
of  the  CVM.  The  pair  exchange  energy  W  for  the  CSA  calculation  was  taken  from  Oates 
et  al.  [990at],  As  shown  in  Fig.  11(b)  the  CSA-calculated  IPB  energies  for  the  (100) 
plane  are  in  accord  with  those  given  in  Fig.  11(a)  calculated  by  Kikuchi  and  Cahn  from 
the  congruent  point  down  to  low  temperatures.  As  also  shown  in  Fig.  12(b),  the  CSA- 
calculated  concentration  profile  for  the  (100)  IPB  between  these  two  phases  at  the 
congruent  point  are  also  in  accord  with  the  CVM-calculated  values  shown  in  Fig.  12(a). 
These  results  led  us  to  believe  that  the  CSA  could  be  used  to  calculate  IPB  energies  and 
should  be  explored  further. 
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Figure  1 1  Calculated  coherent  IPB  energies  of  the  fcc-(Cu,  Au)  and  Ll2-Cu3Au  phases  using  the 
tetrahedron  approximation:  (a)  CVM-Calculation  for  the  (100)  &  (1 10)  planes,  (b)  CSA  calculation 

for  (100)  planes. 


Number  of  the  plane 


Figure  12  Calculated  concentration  profile  for  the  (100)  IPB  at  the  congruent  point:  (a)  CVM  and  (b)  CSA. 


As  noted  previously  [79Kik,  96Astl],  the  vanishing  of  the  calculated  IPB  energy  at  OK 
using  the  tetrahedron  approximation  is  due  to  the  neglect  of  the  longer-range  interaction 
in  the  lattice.  Using  the  tetrahedron-octahedron  approximation  of  the  CVM  (TO-CVM), 
expressed  in  terms  of  the  parameter  a  =  W2/W)  for  a  =  -0.1  and  -1.0,  Asta  [96Astl] 
calculated  the  phase  diagrams  of  the  prototype  Cu-Au  binary  and  then  the  coherent  IPB 
energies  between  the  fcc-(Cu,Au)  and  the  LU-CCujAu)  phases.  The  terms  W2  and  Wi 
denote  the  next  nearest  neighbor  and  the  nearest  neighbor  pair  exchange  energies, 
respectively.  Figures  13(a)  and  (b)  show  the  TO-CVM  and  TP-CSA  calculated  phase 
diagram  and  coherent  IPB  energies  between  the  fcc-(Cu,  Au)  and  LI2-CU3AU  phases  for 
the  (100)  planes.  The  TO-CVM  calculated  values  were  taken  from  Asta  [96Astl]  and 
shown  in  these  figures  as  discrete  points.  As  shown  in  Fig.  13(b)  the  coherent  IPB 
energy  at  OK  is  finite  and  decreases  first  with  increasing  temperature  and  then  change 
slowly  with  further  increases  in  temperatures.  We  used  the  tetrahedron-pair 
approximation  of  the  CSA  (TP-CSA)  to  calculate  the  phase  diagram  and  coherent  IPB 
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Concentration:  X(B)  T/Tc 

Figure  13  (a)  Calculated  prototype  Cu-Au  phase  Figure  13  (b)  Calculated  IPB  energies  between  the 
diagram  with  a  = -0. 1  with  the  TO-CVM  fee  and  LI 2  phases  with  a  =  -0.1.  The  TO-CVM 

calculated  values  shown  as  discrete  points  calculated  values  are  shown  as  discrete  points 

[96Astl  ]  and  the  TP-CSA  values  as  smoothed  [96Astl  ]  and  TP-CSA  values  as  a  smoothed  curve. 


Temperature  (C)  Temperature  (C) 


Figure  14  (a)  Comparisons  of  TP-CSA-calcuiated  Figure  14  (b)  Comparisons  of  TP-CSA-calcuIated 
coherent  IPB  energies  between  fcc-(Al)  and  Ll2-  coherent  IPB  energies  between  y  and  y*  in  Ni-Al 

AI3Li  with  the  TO-CVM  values  and  the  with  experimental  data  [90Mar,  95Ard]. 

experimental  data  [91  Hoy,  84Bau], 

energies  [99Vak,  06Cao],  Our  TP-CSA  calculated  reduced  temperatures  expressed  as  a 
function  of  composition  as  shown  in  Fig.  13(a)  represented  by  smoothed  curves  are  in 
accord  with  the  TO-CVM  calculated  values.  Likewise,  the  TP-CSA  calculated  coherent 
IPB  energies  given  in  Fig.  13(b)  shown  as  a  smoothed  curve  is  in  accord  with  the  TO- 
CVM  calculated  results.  We  did  not  calculate  the  values  down  to  OK  due  to  numerical 
difficulties,  which  will  be  resolved  in  the  near  future.  Nevertheless,  the  fact  that  the 
calculated  values  are  in  reasonable  accord  with  the  TO-CVM  calculated  results  [96Astl] 
is  encouraging. 

We  next  extended  the  TP-CSA  to  calculate  the  coherent  IPB  energies  between  the  fee  and 
LI  2  phases  in  two  real  alloy  systems,  i.e.  Al-A^Li  and  Ni-Al.  As  shown  in  Fig.  14(a)  the 
TP-CSA  calculated  IPB  energies  are  in  accord  with  the  TO-CVM  calculated  values 
[96Ast].  Within  the  uncertainties  of  the  measurements,  there  is  agreement  between  the 
calculated  values  and  the  experimental  data  of  Hoyt  and  Spooner  [91Hoy]  and  Baumann 
and  Williams  [84Bau].  Figure  14(b)  shows  that  the  TP-CSA  calculated  IPB  energies 
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between  y  and  y’  in  Ni-Al  are  in  accord  with  the  experimental  values  of  Ardell  [95Ard]. 
Ardell  re-evaluated  the  data  of  Marsh  and  Chen  [90Mar]  and  obtained  values  in  the  range 
of  2  mJ/m2  but  not  shown  in  this  figure.  It  is  worth  noting  that  Dr.  Chris  Woodward  of 
the  AFML,  Dayton,  OH  [2004Woo]  shared  his  calculation  of  IPB  energies  between  y  and 
y’  in  Ni-Al  with  us.  He  used  a  “first  principles”  cluster  expansion  method  coupled  with 
Monte  Carlo  calculations  to  estimate  the  interfacial  width  and  excess  free  energy  in  y  and 
y’  in  Ni-Al.  For  temperatures  in  the  range  of  experimental  measurements  (see  Figure 
14(b))  his  calculated  IPB  energies  are  in  good  agreement  with  our  estimates,  and  with  the 
available  results  from  experimental  measurements.  The  concentration  profile  across  a 
coherent  (100)  y/y’  interface  obtained  from  recent  grand-canonical  Monte  Carlo 
simulations  at  700K  [04Mis]  is  shown  in  Figure  15(a).  The  concentration  profile 
obtained  from  the  current  study  using  the  TP-CSA  is  also  plotted  in  Figure  15(b)  for 
comparison.  As  pointed  out  by  Mishin  [04Mis],  the  long-range  order  does  not  drop  to 
zero  abruptly  at  the  interface  plane  but  rather  decays  gradually  over  4  -  6  (2  0  0)  atom 
layers  around  the  interface,  as  does  also  the  average  A1  concentration.  The  comparison  in 
Figures  15(a)  and  (b)  shows  that  the  TP-CSA  predicts  a  very  similar  interfacial  width 
with  the  one  derived  from  the  Monte  Carlo  calculations,  although  there  is  a  small 
difference  in  the  average  A1  concentrations  for  the  bulk  y  and  y’  phases,  which  is  due  to 
the  fact  that  the  bulk  phase  diagram  determined  by  these  two  methods  are  slight  different 
from  one  another.  All  these  results  indicate  that  the  CSA,  computationally  less 
demanding  than  the  CVM,  offers  an  alternative  and  perhaps  practical  approach  to 
calculate  (or  estimate)  coherent  IPB  energies  for  y  and  y’  in  other  Ni  binaries  and  even  to 
higher  order  alloys  of  practical  importance. 


Distance  (nm)  Distance  (nm) 


Figure  15  (a)  Concentration  profile  across  the 
coherent  (10  0)  y/y  interface  obtained  by 
grand-canonical  Monte  Carlo  simulations  at 
700  K.  The  points  represent  A1  concentration 
in  individual  (2  0  0)  layers  parallel  to  the 
interface.  The  arrow  marks  the  initial  position 
of  the  interface  at  OK  [04Mis], 


Figure  15  (b)  Concentration  profile 
across  the  coherent  (10  0)  y/y1 
interface  obtained  from  the  present 
study  using  the  TP-CSA. 
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Anti-phase  Boundary  (APB)  Energy  Calculation  -  Cluster/Site  Approximation 


After  extending  the  CSA  to 
multicomponent  systems  and  using  the 
CSA  to  calculate  IPB  energies,  we  use 
the  CSA  to  calculate  anti-phase 
boundary  (APB)  energies  in  Cu-Au 
binary  using  a  similar  approach  as  has 
been  done  for  IPB  energies.  Fig.  16 
shows  the  CSA-calculated  (solid  line) 
and  the  CVM-calculated  (scattered 
points)  {100}  APB  energies  as  a 
function  of  temperature  with  the 
composition  of  Au,  x(Au)  =  0.267. 

Both  results  indicate  an  identical  trend 
and  predict  a  maximum  in  a.  The 
CSA-calculated  value  of  Ga2/KBTc  at 
the  congruent  point  is  around  0.2  for 
the  (100)  APB  and  0.1  for  the  (100) 

IPB.  This  agrees  with  the  relation  between  gApb  and  Oipb  for  (h,  k,  0)  orientations  at  the 
congruent  point  given  by  Kikuchi  and  Cahn  [79Kik],  i.  e.  gApb  =  2gipb.  Thus  the  (100) 
APB  could  either  consist  of  two  (100)  IPB’s  with  a  disordered  layer  in  between  or  it 
coincidently  has  the  same  free  energy.  Fig.  17(a)  and  17(b)  show  the  calculated  profiles 
of  the  APB  and  IPB  and  conform  that  the  disordered  layer  is  part  of  the  APB  and  that 
indeed  the  APB  has  become  two  IPB’s.  As  the  temperature  decreases  from  the  congruent 
point  to  T/Tc  =  0.9,  the  APB  narrows  but  the  internal  structure  of  the  APB  remains  the 
same  as  shown  in  Fig.  17(a).  However,  according  to  our  CSA-calculated  results,  there  is 
a  structure  change  when  temperature  decreases  to  T/Tc  =  0.83.  This  is  also  in  accord 
with  those  results  obtained  from  the  CVM  [79Kik].  But  the  CSA  calculation  is 
computational  much  less  demanding. 


Figure  16.  Comparison  of  (100)  APB  energies  as  a 
function  of  temperature  at  the  constant  composition 
x(Au)  =  0.267  calculated  by  CVM  and  CSA 


Figure  17.  The  calculated  profiles  of  the  (100)  APB  (a)  and  IPB  (b)  at  the  congruent  point. 
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Number  of  the  plane  Number  of  the  plane 

Figure  1 8.  There  is  a  structural  change  in  the  core  of  (100)  APB’s  between  T  /  Tc  =  0.9  (a)  and  0.83  (b). 
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Chemistry  and  Materials  Properties/Oxidation  and  Corrosion  at  the  Gordon 
Conference  on  High  Temperature  Materials,  Processes,  and  Diagnostics  2006, 
Colby  college,  ME,  USA,  07/19/06.  (“Application  of  the  CSA  to  FCC  Phases  in 
Calculating  Multi-Component  Nickel-Based  Alloys”,  Presented  by  W.  Cao). 
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The  following  four  have  been  awarded  to  Y.  Austin  Chang 


2)  The  2005  Distinguished  Award,  International  Workshop  on  Advanced 
Intermetallic  and  Metallic  Materials. 

3)  Wisconsin  Idea  Fellow,  University  of  Wisconsin  System,  2004-2005. 

4)  Distinguished  Speaker,  Department  of  Materials  Science  &  Engineering,  Penn 
State  University,  2004. 

5)  Identified  as  a  highly  cited  materials  scientists  by  ISIHighlyCited,  2003,  covering 
the  period  1981-1999. 

6)  Edward  DeMille  Campbell  Memorial  Lecturer  Award,  ASM-Intemational,  2003. 
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